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A brief review is given of the sign problem in finite density lattice QCD and various attempts to overcome it. To 
date there is still no solution to this problem which would work for realistic QCD. The main focus then is on the 
deconfined phase, where QCD can be described by a dimensionally reduced effective action. After summarizing 
derivation and validity of the effective theory, it is demonstrated that it can be simulated efficiently in the presence 
of a chemical potential for quarks fi/T < 4. Direct comparison of simulations with imaginary and real /i suggests 
that equilibrium plasma properties could be analytically continued from 4d QCD simulations at imaginary [i. 



1. INTRODUCTION 

QCD at finite baryon density occurs in nature 
in two rather different regimes: i) In heavy ion 
collisions, whose initial state has non-zero baryon 
number, and in which any subsequent plasma 
state is a state of high temperature and low den- 
sity, ii) In the core of neutron stars, which is 
composed of cold and very dense matter. These 
two situations correspond to the regions close to 
the axes in the tentative QCD phase diagram 
Fig. [j]. While in the latter case some interesting 
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Figure 1. Conjectured phase diagram for QCD 
with two massless and one heavy flavour. From 
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conjectures about a colour-superconducting state 
of matter have been made understanding the 
former regime is particularly pressing in view of 
running and upcoming experiments at SPS, LHC 
(CERN) and RHIC (BNL). Eventually we would 
of course like to know the structure of the en- 
tire phase diagram, in particular the position and 
nature of the transition lines separating different 
regions and the (tri-) critical points. 

In order to come to first principles predictions, 
non-perturbative methods are required. Unfortu- 
nately, lattice QCD has so far failed to be a vi- 
able tool for the analysis because of the so-called 
"sign-problem". The QCD partition function is 
given by 



Z 



DU det M (p)e 



-s g [u] 



S f = j d 4 x 

colour, flavour 

M = £>-^7o, V = j-(d + iA)+m, (1) 
where S g is the pure gauge action. The relation 

75 Af 75 = Aft (2) 

implies positivity of the fermion determinant for 
fi = 0, detM(0) > 0. However, for SU(3) and 
/i 7^ 0, Eq. (||) does not hold, detM is complex 
and prohibits standard Monte Carlo importance 
sampling, for which a manifestly positive measure 
in the functional integral is required. 

In the first part of this contribution, a short 
review is given on the various attempts and tech- 
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niques proposed in the literature to deal with this 
problem. Many of these try to implement a chem- 
ical potential at zero or low temperature and to 
obtain the "onset" value fj, = /^o, where the transi- 
tion from a zero density state of matter to a state 
of finite density takes place. At zero tempera- 
ture this should happen when one baryon can be 
made, i.e. [io = tob/3. Since most of these tech- 
niques have been reviewed in detail in previous 
years [HQl, this section will be rather brief. 

The larger part is then devoted to the study 
of the plasma situation at high temperature and 
low density, where it was recently demonstrated 
that numerical simulations of static quantities at 
finite density are feasible in the framework of a 
perturbatively derived effective theory g. 

2. NUMERICAL TECHNIQUES FOR 
FINITE DENSITY 

2.1. Reweighting methods 

The Glasgow method |j| evaluates the partition 
function by absorbing the complex determinant 
into the observable, and doing the importance 
sampling with the positive part of the measure, 



(O) 



OdetM(fi) 
detM(// = 0)/„ = 



(3) 



One may further expand the grand canonical par- 
tition function in the fugacity z = exp fi/T, 



Z(pi,T, V) = 5> B Z(i?,T,^), 



(4) 



where Zb is the canonical partition function at 
fixed baryon density, and compute the terms in- 
dividually. However, numerical results give for 
the onset of baryon density the unphysical value 
Ho = instead of the expected (Xq — 

which is also the pathological behaviour of the 
quenched theory ||. It seems now understood 
that the reason for this failure is the poor overlap 
between the ensemble of configurations at // = 
used to compute averages by reweighting, and the 
correct ensemble at [i ^ 0. 

This situation can be improved by splitting 
the determinant in modulus and phase, det M = 
| det Af | exp(i(/)), and to incorporate the modulus 



into the measure, while only the phase is used for 
reweighting Q, 
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The expectation value of the phase factor can be 
viewed as a ratio of two partition functions, one 
of the true theory and one with the positive mod- 
ulus of the determinant in its measure. It behaves 
thus as exp — AF, where the exponent is the free 
energy difference and an extensive quantity. Con- 
sequently, (exp i(f>) ~ exp(V^), which requires ex- 
ponentially large statistics for realistic volumes. 
A similar variant of the Glasgow method works 
for heavy quarks g. 

2.2. Static quarks 

The sign problem can be handled in a static 
quark model, in which all quarks are infinitely 
heavy ||[l0) ■ It is obtained by sending both quark 
mass and chemical potential to infinity /i, m — > 
oo, while fixing a finite baryon number. This is 
achieved by keeping the constant 

C = (2am/exp{a^)) Nt (6) 

fixed, where N t is the temporal extent of the lat- 
tice. In this case an analytic expression for the 
determinant is obtained which is trivial to evalu- 
ate. The available computing time then suffices to 
gather large enough statistics for the phase of the 
determinant to be calculated with sufficient accu- 
racy. The question is whether results obtained for 
static quarks translate at all to QCD with small 
quark masses. A more detailed discussion of this 
approach is presented in another contribution to 
this conference 

One may attempt to extend this approach to 
finite but large quark masses by going beyond 
the leading term in the 1/m-expansion, based 
upon an interesting observation ]l2| |: analysis of 
the data of [jlO| shows a correlation between the 
phase of the fermion determinant and the imagi- 
nary part of the Polyakov line, cf> ~ VPi, as shown 
in Fig. ^| Indeed, the leading order term in 1/m 
for the phase is (ll| 

cj) = 2(2/m) Nt Vsmh(N t n)Pi . (7) 

The correlation is weakened by higher terms in 
the quark mass expansion which do not show 
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Figure 2. Correlation between the phase of the 
fermion determinant and the imaginary part ot 
the Polyakov line for static quarks at finite den- 
sity on various volumes. From [p^[ . 



this proportionality, and therefore it weakens at 
smaller quark masses. 

This suggests the following approximation: 
computing a constrained partition function 

Z' = J DU det Me~ s ^5(Pi[U}) (8) 

eliminates configurations with significantly fluc- 
tuating phases <f> whenever its correlation with 
Pi is strong, and thus improves the sign prob- 
lem significantly. The approximation is good, 
i.e. Z' « Z, if Pi -c P r , which is the case at 
small n/T, or if P r w P t « at T - 0, fj, < /Lt . 
We shall come back to this point when hot QCD 
at finite density is discussed. 

2.3. Imaginary chemical potential 

There are a number of suggestions to consider 
imaginary chemical potential, for which the inte- 
gration measure is positive and simulations are 
straightforward. The relation to real chemical 
potential is provided by the canonical partition 
function at fixed quark number Q, which is re- 
lated to the grand canonical partition function at 
fixed imaginary chemical potential by Jl3[ ] 

Z(T, Q) = C T % Z(T, /i = i M/ ) e-*'Q' T .(?) 
J-ttT 2nl 



One proposal is to simulate Z for various values of 
^i, and then numerically do the Fourier integra- 
tion This becomes more and more difficult, 
however, for large Q, and extrapolation to the 
thermodynamic limit seems questionable. The 
method has been tested in the two-dimensional 
Hubbard model @, but not for QCD. 

A second proposal is to compute expectation 
values with the grand canonical partition func- 
tion Z(i/i,T), fit them by power series in [i and 
analytically continue the series to real fi. This 
has been explored for the chiral condensate in the 
strong coupling limit |L5| . 

Finally, starting from Eq. (Jjj), a well-defined 
quenched limit of the theory at finite density can 
be obtained Jig . The limit m — * oo is taken in 
the expression Eq. (||), upon which it goes over 
into 

Z(T,Q)= J DU f Q [P] e - s °M, (10) 

where /q[P] is a functional of Polyakov lines. 
The quenched limit thus is not simply realized 
by replacing the fermion determinant by a con- 
stant. Instead, there are remainders of the static 
quark propagators in the partition function, the 
Polyakov loops, which carry the baryon number. 
The expression after the limit is still complex, 
and thus continues to have a sign problem. How- 
ever, simulations have shown that it is weakened 
enough to be well manageable numerically, and 
results for the static potential are available |l6f| . 

2.4. Avoiding the sign problem: SU(2) 

In view of the difficulties with SU(3), one may 
resort to studying two colour QCD which, thanks 
to the realness of the SU(2) matter representa- 
tions, has a real and positive fermion determi- 
nant and can be simulated straightforwardly []l7f . 
Unfortunately, the same property leads to other 
peculiarities which make it less analogous to QCD 
than one might have hoped: charge conjugation is 
related to global SU(2) rotations. Thus, the con- 
densates (qq) and (qq) are related by rotations 
as well, and there is no fundamental distinction 
between them. Since there are only two colours, 
'baryons' consist of even numbers of quarks, and 
through SU(2) rotations they are equivalent to 
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mesons. In particular, this results in a physi- 
cal onset chemical potential fio = /2, even in 
the unquenched theory. Moreover, since (qq) is 
gauge-invariant, there can be no colour supercon- 
ductivity either. Hence, the model is simulable, 
but some main questions of interest for QCD have 
disappeared. 

Somewhat surprisingly, these questions reenter 
if the representation of the fermions is taken to be 
adjoint and only one flavour is considered pq . In 
this case, there are no qq Goldstone modes. More- 
over, the theory permits gauge-invariant three- 
quark baryon operators qqq with spin 1/2, and 
thus no for the finite density onset should be dif- 
ferent from half the meson mass. Further, (qq) is 
gauge- variant and hence permits a study of colour 
superconductivity. Finally, in the lattice model 
with an adjoint fermion, the determinant is real 
but not positive definite, hence a weakend sign 
problem is present which can be handled by a 
two-step multi-boson algorithm fjg|] . 

More detailed discussions of SU(2) models to- 
gether with numerical results may be found in the 
finite density section of these proceedings . 

2.5. Finite isospin density 

In has also been proposed to consider a chem- 
ical potential for isospin rather than for baryon 
number |l4"|]2C| l. In this case, and with two de- 
generate light flavours m u = md, Eq. (||) can be 
supplemented by an isospin rotation to give 



Ti7 5 Af7 5 Ti = Af f , 



(11) 



which again implies det M > and thus can 
be simulated. Analytic work for very small and 
large isospin chemical potential, where chiral per- 
turbation theory or ordinary perturbation theory 
are applicable, respectively, lead to the conjec- 
ture that the hadron and quark matter phases 



are smoothly connected 20 



2.6. Hamiltonian formulation 

Another attempt addresses the problem in a 
Hamiltonian formulation of lattice QCD at fi- 
nite density |2|]. Quark number density, chiral 
condensate, some meson and the nucleon masses 
have been calculated in the limits of a free and a 
strongly coupled theory at zer temperature. The 



chiral phase transition is found to be first order, 
and the onset fi$ is consistent with ma/3. So far, 
however, no numerical implementation of this ap- 
proach is available. 

2.7. Cluster algorithm 

Finally, a new algorithm has been proposed to 
actually solve the sign problem p2]| . So far, it 
has only been worked out for non-relativistic spin 
systems. The approach uses cluster algorithms 
to decompose a fermion configuration into clus- 
ters of spins that can be flipped independently, 
individually changing sign or not under the flip. 
This enables to cancel negative against positive 
contributions, leaving only net positive contribu- 
tions. Details and possible extensions to QCD are 
discussed in the contribution |ll| . 

2.8. Finite temperature and density 

All proposals that have been implemented nu- 
merically have thus failed on realistic QCD so far, 
and the failure is particularly drastic at low tem- 
peratures. The methods that do work numeri- 
cally are limited to static or very heavy fermions, 
and they work better at high temperatures, sug- 
gesting that the sign problem may be easier to 
overcome in hot QCD. Let us then shift attention 
to temperatures above the deconfinement transi- 
tion. 

3. EFFECTIVE THEORY APPROACH 
TO HOT QCD 

For temperatures larger than a few times the 
deconfinement temperature T c , the static (equi- 
librium) physics of QCD can be described by an 
effective theory, which indeed permits simulations 
with non- vanishing real chemical potential |j] ■ Of 
course, the restriction to the deconfined phase 
presents a limitation. Another one is that finite 
density simulations are possible only up to some 
value of the chemical potential, where the sign 
problem sets in. Presently the accessible range 
is /i < AT. Finally, if the effective theory is to be 
quantitatively accurate, one has to be sufficiently 
away from the deconfinement transition. 

If there are limitations, there are also strengths 
of this approach. The effective theory is three- 
dimensional and purely bosonic, so simulations 
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are cheap and yield rather precise continuum lim- 
its, leaving little ambiguity in the interpretation 
of the results. The approach is good for SU(2) 
as well as SU(3), and permits to study any num- 
ber of quarks with any small or zero mass. Fi- 
nally, the domain of parameter space where the 
approach is applicable contains the phenomeno- 
logically relevant region in which heavy ion colli- 
sions are operating. At and above SPS energies, 
densities in heavy ion collisions are estimated to 
be Hb/T < 4.0 p3[ ], i.e. a quark chemical potential 
yit/T< 1.3, which is well within the range where 
simulations are feasible. 

3.1. Equilibrium correlation lengths 

Let us first discuss the physics question that 
will be addressed in the following. Consider 
gauge- invariant, local operators A(x). Static 
equilibrium physics is described by euclidean time 
averages 



l/T 

A(x) = T I drA(x,-ir). 



(12) 



The spatial correlation function of such operators, 

C(M) = (A( x )A(0)) c ~e- M M, (13) 

falls off exponentially with distance. The "screen- 
ing masses" M have a precise non-perturbative 
definition: they are the eigenvalues of the space- 
wise transfer matrix in the corresponding lat- 
tice field theory. Physically, they correspond to 
the inverse length scale over which the equili- 
brated medium is sensitive to the insertion of a 
static source carrying the quantum numbers of 
A. Beyond 1/M, the source is screened and the 
plasma appears undisturbed. With these def- 
initions, all screening lengths corresponding to 
gauge-invariant sources can be computed on the 
lattice in principle. 

3.2. Dimensional reduction 

If we are interested in length scales larger than 
the inverse temperature, |x| ~ 1/gT ^s> l/T, the 
situation simplifies considerably. In this case the 
integration range for euclidean time averaging in 
Eq. ( |l2| ) becomes very small, and the problem 
effectively three-dimensional (3d). 

The calculation of the correlation function 
C(|x|) can be factorized: the time averaging 



may be performed perturbatively by expanding 
in powers of the ratio of scales gT/T ~ g, which 
amounts to integrating out all modes with mo- 
menta ~ T and larger, i.e. all non-zero Matsubara 
modes, in particular the fermions. This procedure 
is known as dimensional reduction J24J]. It is in 
the spirit of a Wilsonian renormalization group 
approach, where an effective action for coarser 
scales is derived by averaging over the smaller 
scales. However, here this step is done perturba- 
tively and entails two approximations: the first 
is that the computation is obviously limited to a 
finite order in g. The second is to neglect higher- 
dimensional operators, which are suppressed by 
powers of the scale ratio^J. 

The remaining correlation function of 3d fields 
is then to be evaluated with a 3d purely bosonic 
effective action, which describes the physics of the 
modes ~ gT and softer. The screening masses 
M now correspond to the eigenvalues of the 3d 
transfer matrix. Simulation of this reduced prob- 
lem is easy. Without fermions and one dimension 
less, much larger volumes and finer lattices can 
be considered. Moreover, 3d gauge theories are 
superrenormalizable and the coupling scales lin- 
early with lattice spacing. Hence, very accurate 
continuum limits can be obtained. 

The effective theory emerging from hot QCD 
by dimensional reduction is the SU(3) adjoint 
Higgs model p4p^p7|| with the action 



S = 



+mrrr(^) + A 3 (Tr(^) 2 }. (14) 



As 4d euclidean time has been integrated over, 
Aq now appears as a scalar in the adjoint repre- 
sentation. The associated Hamiltonian respects 
SO (2) planar rotations, two-dimensional parity 
P, charge conjugation C and A -reflections R. 

The parameters of the effective theory are via 
dimensional reduction functions of the 4d gauge 
coupling g 2 , the number of colours N and flavours 
Nf, the fermion masses and the temperature T. 
In all of the following fermion masses are assumed 



2 A method to construct an effective theory by a (pertur- 
bat ive) block spin transformation has been proposed in 
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to be zero, but in principle any other values may 
be considered as well. At leading order in the 
reduction step one has 



9l 



m 3 
A 3 



N N f \ 

(6 + N — Nf)g i (fl)T . 



1/2 



24tt 2 



(15) 



Specifying T/A m , N f 



Note that is just the leading order Debye mass 
pq| . The reduction step has been performed up 
to two-loop order p'jj ] at which parameters have 
relative accuracy 0(g 4 ) and 1-loop scale depen- 
dence on jl is cancelled 
then completely fixes the effective theory. 

If we are interested in the infrared-sensitive, 
non-perturbative scales of order ~ l/g 2 T p9| , 
the fields A ~ gT may be integrated out as well, 
leaving us with a 3d pure gauge theory. In order 
to decide which is the correct effective theory for 
a given temperature beyond parametric consid- 
erations, simulation results for screening lengths 
will be compared between the reduced and the 
full theory. 

With the simulation part being rather accurate, 
the main error in the correlation functions is due 
to the reduction step, and at two-loop level this 
is formally of the order |BE 



(16) 



In the treatment of the electroweak phase transi- 
tion, this is much less than 5% [|o), but for hot 
QCD the coupling and thus the error is larger. 
In fact, one may wonder why perturbative reduc- 
tion should work at all, given that the effective 
g 2 ~ 2.7 at T = 2T C @. In the loop expan- 
sion, the coupling appears as a s /ir rather than g. 
Inspecting the computed series for the effective 
theory parameters Eq. (15) up to next-to- leading 
order, it appears that it converges surprisingly 
well at temperatures not much above the critical 
one |27 32 1. The other source of error is neglect- 
ing higher dimension operators. This is justified 
as long as the dynamical mass scales of the ef- 
fective theory are smaller than those integrated 
out ~ 2ttT. The success of the approach then 



15 



N f =0, T=2T C , L=40, (3=28 



t; 



X 


5 5 


I ■ I O 


■ 5 






o 







oy oy oy o _ 1 + 1_ 2 

N,=0, T~10"T C , L=40, p=28 



1.5 



1.0 



7 i ? 



i 

0.5 

o.o L-^—^—^ 



I t 



1 1 



Figure 3. Spectrum of screening masses in var- 
ious quantum number channels at low and high 
temperatures, for SU(3), Nf = 0. Filled symbols 
denote 3d glueball states. From 0. 



depends on whether there is a non-perturbative 
scale separation. This can only be checked after 
simulating the effective theory. 

3.3. Numerical results 

Results from simulations || of the gauge- 
invariant screening spectrum with the effective 
theory Eq. ( [lij) are displayed in Fig. [|. Screening 
masses are classified according to the symmetries 
of the 3d Hamiltonian by the quantum numbers 
Jr C ■ The open symbols show states dominated 
by scalar operators such as Tr(Ao), Tr (AoFf 2 ) 
etc., whereas full symbols show states receiving 
only gluonic contributions Tr(F : 2 2 ) etc. It is a re- 
markable finding of detailed mixing analyses in 
several models [p3]]34|,|5| , that the latter are quan- 
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titatively consistent with the 3d glueball states 
|35[ and apparently completely insensitive to the 
presence of the Aq scalar field. Because of this 
pronounced non-mixing, we can thus state a first 
important result: For any reasonable temper- 
atures the largest correlation length of gauge- 
invariant operators belongs to the Ao degrees of 
freedom and not to the Ai , in contrast to the naive 
parametric picture. This demonstrates that the 
physics from the scale gT down to g 2 T is com- 
pletely non-perturbative, and hence A may not 
be integrated out perturbatively. Only for tem- 
peratures T > 10 7 T C becomes the coupling small 
enough that the parametrically expected ordering 
of screening masses is realized. 

3.4. Testing dim. red. in hot Yang-Mills 
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Figure 4. Comparison of screening masses 
(C — +) in hot SU(2) pure gauge theory as deter- 
mined in 4d (empty symbols) and 3d (full sym- 
bols) effective theory simulations. From J34|. 



We now wish to know the accuracy of the re- 
duced theory. Fig. || compares the results for hot 
SU(2) gauge theory as obtained in the 4d lat- 
tice theory |3(| with those from the effective the- 
ory |Q. Note that the effective theory is only 
valid up to its cut-off M/T ~ 2ir and above this 
level disagreement is to be expected. The appar- 
ent excellent agreement for the two lowest lying 
states should be taken with some care given some 



Table 1 

Comparison of the lowest two screening masses in 
SU(3) Yang-Mills at T = 2T C , obtained in 3d § 



and 4d 



56l simulations. 



M/T 


3d 


4d 


0++ 


2.71(6) 


2.60(4) 


01+ 


4.6(2) 


6.3(2), 



systematic uncertainties. In 4d simulations, the 
projection properties of operators and the scal- 
ing behaviour are more difficult to control than 
in 3d. In the effective theory, there is an ambi- 
guity in the choice of renormalization scale. The 
direct comparison may thus still suffer from ef- 
fects up to 20%. For the case of SU(3) one finds 
again quantitative agreement in the largest corre- 
lation length but indeed about 20% deviation in 
the next shorter one, cf. Table § Good agreement 
between the 4d and the 3d effective theory is also 
found with other observables, like the static po- 
tential |^6| , the spatial string tension j3?j] , the De- 
bye mass Q and gauge- fixed propagators |3^] .[] 
Thus, dimensionally reduced pure gauge theory 
gives at least a semi-quantitative description of 
the largest correlation lengths in the system down 
to temperatures as low as ~ 2T C . 

3.5. Inclusion of fermions 

The main advantage of the effective theory is 
that fermions, having always non-zero Matsubara 
frequency, are treated analytically. Any change in 
the number of fermion species or their masses is 
encoded in the parameters of the effective theory, 
Eq. (|l5|). The results for the spin zero channels 
at various Nf (all massless) are shown in Fig. ||. 
The screening masses get larger with increasing 
Nf , so that even the smallest of them comes close 
to the cut-off ~ 2ttT of the effective theory, thus 
diminishing the scale separation needed for the 
effective theory to be valid. Consequently, only 
the lowest state may be expected to describe hot 
4d physics. As in the Nf = case, this situation 
is improved by going to larger temperatures. 

Another, and more severe, limitation comes 



3 Dimensional reduction also works to high precision from 
3d finite T to 2d |40| . 
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Figure 5. Spectrum of scalar screening masses 
for various Nf with zero fermion masses. Filled 
symbols denote 3d glueball states. From ||. 
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Figure 6. The smallest scalar screening masses 
due to correlations of gluonic (g) and mesonic (it) 
operators. The lines represent phenomenological 
fits. From ElL 



from the fact that, close to T c , the fcrmionic 
modes begin to feel the chiral phase transi- 
tion and thus become very light through non- 
perturbative effects. At some point this effect 
will be so large that they constitute the light- 
est degrees of freedom and may no longer be in- 
tegrated out. This situation has been demon- 
strated recently in a 4d simulation with Nf = 4 
light fermions jyj , and is shown in Fig. |[ Below 
T ~ 2 — 3T C , the pion correlator yields the longest 
correlation length, and the purely bosonic effec- 
tive theory fails to give the correct infrared be- 
haviour. This problem becomes worse for larger 
Nf, whereas it is slightly better for the realistic 



case of two light fermion flavours. It would be 
very interesting to see if this behaviour can be re- 
produced in the effective theory, if NRQCD meth- 
ods are applied and fermionic degrees of freedom 
are kept in the effective action |ft^| . 

4. DIMENSIONAL REDUCTION AT FI- 
NITE DENSITY 

Having seen that dimensionally reduced QCD 
gives a reasonable description for the longest cor- 
relation length at temperatures T>2 — 3T C , we 
may now return to the problem of introducing a 
chemical potential ||. The contribution of the 
fermion determinant with [i ^ to the effective 
potential of A has been evaluated at the two- 
loop level j43|. From this one may derive the 
corrsponding changes in the three-dimensional ef- 
fective action, brought about by the chemical po- 
tential term. To leading order, it generates one 
extra term in the action Eq. (jl4|) and changes the 
bare Aq mass parameter, 

(PxTiAq, z = — gT^J 

( M \ 2 3N f 
\ttTJ 2N + N f _ 



S 



m>3 



S- 



m>3 



1 



(17) 



The extra term S z is odd under R, C, and hence 
the action no longer respects these symmetries, 
while parity is left intact. Consequently, screen- 
ing states now are labelled by J p . 

The effective action is complex just as the 4d 
original one, and hence there is a sign prob- 
lem. Expectation values have to be computed by 
reweighting with the complex piece of the action 

-ixS x \ 



(O) 



(Oe 



(18) 



(cos(z5^)) ' 

where the subscript denotes averaging with /i = 0. 
Cancelling contributions to the expectation value 
occur whenever (zS z ) 3> 1. The perturbative 
width of the distribution at two-loop level can 
be computed from the Aq sunset diagram to give 

. , _ , z /6L\3/2 



In — + 0.58849 ) 

lzarriD J 



1/2 



(19) 
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Figure 7. Distribution of the reweighting factor 
{zS z ) for fixed /i/T and various L (top), and for 
fixed L and various fi/T (bottom). From ||. 
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Figure 8. The lowest scalar screening masses 
when an imaginary (top) or real (bottom) chem- 
ical potential is switched on. From ||. 



One recognizes the same behaviour as in the 4d 
case, in that the width of the phase distribution 
is growing with ~ VV(fi/T). Fig. ^ shows some 
distributions of (zS z ) as obtained by Monte Carlo 
for various volumes and chemical potentials. As 
long as /1 < 4T, the distribution is well contained 
within [— 7r, it) for volumes large enough so that 
the masses attain their infinite volume values || . 
Correspondingly, as long as this condition is met, 
there are no problems with cancelling contribu- 
tions to correlation functions, and the signal is as 
good as at [i — 0. For even larger values of the 
chemical potential the sign problem sets in and 
statistical errors explode. 

Next, one has to check the overlap between the 



ensembles at \i ^ and /j, = in order to avoid 
a wrong biasing of the importance sampling by 
the reweighting procedure. For a given \/j,\/T, one 
can get an idea about this using imaginary chemi- 
cal potential, where reweighting can be compared 
with the full action update. This is shown in the 
top panel of Fig. |[ Within the small statistical 
errors, no difference between reweighted and orig- 
inal ensembles is noticable. This result then gives 
us confidence in the corresponding ones obtained 
for real chemical potential, bottom of Fig. [| 

Because of the different qualitative behaviours 
of 3d gluonic and scalar states mentioned earlier, 
we may expect to observe a change in the na- 
ture of the ground state excitation at some real 
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value of fi. Indeed, Fig. || bottom suggests a level 
crossing at n/T ~ 4.0. This would mean that 
the longest correlation length in the thermal sys- 
tem does not get arbitrarily short with increasing 
density, but rather stays at a constant level. Note 
that the value of M /T at this crossing is already 
so large that the effective theory may be inac- 
curate quantitatively. However, the qualitative 
effect should be the same in the 4d theory. 

5. RELATION TO 4D 



One may now ask why the sign problem in the 
effective theory is weakened, and what the rela- 
tion to the original problem in four dimensions 
is. To this end, note that the fermionic Matsub- 
ara modes at finite temperature have propagators 
(Piji + m + 7o(2n + l)7rT) _1 . For large T these 
modes are heavy, which is why dimensional re- 
duction integrates them out. We can then make 
contact to the heavy quark expansion of Sec. 2.2, 
where a correlation between the phase of the 
fermion determinant and the imaginary part of 
the Polyakov loop has been observed (cf . Fig. ||) . 
Expanding the exponential of the imaginary part 



channel 



of the Polyakov loop Pi ~ ImTr( J" drA ) 3 + . . . in- 
dicates that after integrating over euclidean time 
the leading contribution is 



T' 



(20) 



which is just the extra term in the effective action 
Eq. ([l?]). There, part of the path integral has 
been computed already analytically, namely the 
integration over the hard modes of order p > T. 
This has produced the numerical factor in z, 
Eq. ( |I^ , |l9| ), suppressing the width of the distri- 
bution relative to the one before this integration. 

6. IMAGINARY VS. REAL CHEMICAL 
POTENTIAL 

With an effective theory permitting simulation 
of real and imaginary chemical potential at hand, 
one may now return to an earlier suggestion and 
study the feasibility of analytic continuation of 
observables from imaginary to real /i Q]. Away 
from phase transitions, the screening masses have 
to be analytic in fi/T, as there are no massless 
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Figure 9. Fits of Eq. fl21| ) to the lowest two 
screening masses at imaginary and real chemical 
potential. From Ej. 



modes in the theory. Moreover, since a change 
fi — * —/j, can be compensated for by a field re- 
definition Ao — * — Aq in Eq. (p"4|), all physical 
observables must be even under this operation. 
In the original 4d theory the same statement fol- 
lows from compensating fi —fj, by a C (or CP) 
operation. For small values of n/T, the screening 
masses may thus be written as 



M 
~T 



Co + Cl 



JL 

nT 



-c-i 



JL 

ttT 



-°(&) -< 21 » 



One may now ask to what extent a truncated 
power series like this reproduces the data. In 
the range where this is possible for both real and 
imaginary fj,, the possibility of analytic continua- 
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tion is easily checked by examining whether the Ci 
agree between the two data sets. Fig. ^| shows the 
lowest states in the system for real and imaginary 
potential, together with the fits to Eq. (21). 

Note that in the imaginary \i case, data are 
only available up to \n\/T ~ 1.5, and hence C2 is 
not well constrained. The reason is that imagi- 
nary chemical potentials favour a Zjv-broken min- 
imum over the symmetric one in the 4d effective 
A potential, once \/j,\/T > 7r/3. Thus a phase 
transition occurs and analyticity is lost. In the 
effective theory, this transition shows up by a de- 
crease of the scalar bare mass parameter sending 
the system into a broken phase. In the region up 
to this critical value, however, complete compat- 
ibility for the values of co,i is observed between 
the two data sets Q . 

One may object that, once the deconfinement 
transition line is continued to complex values of 
/x, it may get distorted, and a continuation might 
fail if it pierces the transition sheet T c (fi) while 
rotating from imaginary to real values. In fact, 
for imaginary /i screening masses are decreasing. 
Correspondingly we would expect the line T c (fi) 
in the {T, Im(/i)}-plane to be an increasing func- 
tion of , which is the distortion we were worried 
about. However, the change in T c ([i) should be 
of the same order in /x/ (irT) as the change in the 
screening masses. For the maximal \n\/T < 1.3 
allowed, the change should not be more than 20%, 
and the procedure is surely safe for T > 1.5T C . 

In the region where the effective theory is appli- 
cable, it is thus found that direct analytic contin- 
uation does seem to provide a working tool for 
determining correlation lengths. This suggests 
that in 4d simulations analytic continuation from 
imaginary chemical potential would give physical 
results if a good ansatz for the tt-dependence is 
available, and would allow to go closer to T c de- 
termining, e.g., the free energy density and the 
spatial correlation lengths there. 

7. CONCLUSIONS 

In summary, numerical simulations of SU(3) 
QCD at finite density are feasible only for static 
or heavy quarks, where the sign problem persists 
but enough computational power is available to 



calculate the phase of the fermion determinant 
with sufficient accuracy. However, all numerical 
approaches so far fail for realistic QCD. 

Qualitative insight may be gained from simu- 
lations of two-colour QCD, where the model with 
one flavour of adjoint fermions shares several fea- 
tures with QCD, or from simulations of QCD at 
finite isopsin density, in which cases there is no 
sign problem. It will be interesting to see fur- 
ther developments of another approach, which at- 
tempts to actually solve the sign problem by use 
of cluster algorithms. 

At temperatures above the deconfinement 
transition, the equilibrium physics of QCD is 
described by a perturbatively derived, three- 
dimensional and purely bosonic effective theory, 
which is the more accurate the higher the tem- 
perature. Reasonable results for the longest cor- 
relation length may be obtained for temperatures 
>2-3T c . 

The effective theory can be efficiently simulated 
even after introduction of a finite chemical poten- 
tial for light or massless quarks, as long as \i < 4T. 
This domain of applicability covers the parameter 
values which are phenomenologically relevant for 
heavy ion collision experiments. It is furthermore 
found that screening masses can be computed by 
analytic continuation of calculations at imaginary 
chemical potential. This points to the possibility 
that static plasma properties could be computed 
in four-dimensional QCD simulations employing 
imaginary chemical potential, thus circumventing 
the sign problem. Another problem for future 
work might be to improve the effective theory by 
NRQCD methods, allowing also hadronic observ- 
ables to be computed. 
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